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Abstract 

A Multi-Scale representation of the combustor dynamics was used to create a self-tuning, scalable controller to 
suppress multiple instability modes in a liquid-fueled aero engine-derived combustor operating at engine-like 
conditions. Its self-tuning features designed to handle the uncertainties in the combustor dynamics and time-delays 
are essential for control performance and robustness. The controller was implemented to modulate a high-frequency 
fuel valve with feedback from dynamic pressure sensors. This scaleable algorithm suppressed pressure oscillations 
of different instability modes by as much as 90 percent without the “peak-splitting” effect. The self-tuning logic 
guided the adjustment of controller parameters and converged quickly toward phase-lock for optimal suppression of 
the instabilities. The forced-response characteristics of the control model compare well with those of the test rig on 
both the frequency-domain and the time-domain. 
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Nomenclature 

2 m><2 m “state-space” matrix in the finite-difference equation of m acoustics modes (Eq. (la)) 
control loop-time, also taken as time step in the “state-space” equation of the “plant” model 
discrete time index associated with T 
valve-commands for fuel modulation (in Volts) 

expected valve effectiveness on fuel-flow modulation (in 1/Volt, so that “Zz is dimensionless) 
index of the modeled modes, starting with “z =1” for the fundamental mode 
approximate natural frequency (in radians/second) of combustor acoustics mode i 
approximate damping coefficient for acoustics mode i (“critically damped” is when “C= 1”) 

2 m-vector of m pairs of “quadratures” (dimensionless, modal state-variables) for m acoustics modes 
2 m-vector function of many variables to model heat-fluctuation lumped effects on acoustics 
(Eq. (la)) 

dimensionless gain parameter for the effect of fuel-flow perturbations on heat-release 
time-shift operator, also interchangeably denotes the associated (possibly random) delay-time 
wide-band, random disturbances 

reference pressure differential, in psi, comparable to the pressure drop across the injector 
combustion zone pressure oscillations (about a mean), in “psi” (pounds per square inch) 

2 m-vector of m pairs of “quadratures” extracted from feedback signals {£„} using quadratures filters 
2 m-vector of state-variables in the dynamic equation to project X n ahead N time-steps, given \S‘ n } 
a prediction for X„ by a Kalman estimate of q n using S l „ as measurement feedback 


I. Introduction 

Lean-burning front-ends (e.g., lean-burn direct injection) for combustors offer the possibility to drastically 
reduce particulate and NO x emissions by aero-engines. The dynamics of this type of combustor, however, are 
susceptible to instability as has been observed in ground-based lean-burning power generators. This is due to the 
reduced internal damping and the increased perturbations from more energetic mixing of liquid-fuel and air. The 
more homogeneous acoustic medium in a typical lean-burn combustor is favorable for well-formed acoustics modes. 
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In addition, reduction or elimination of film-cooling liners in the combustor will result in acoustically hard boundary 
conditions which reflect the perturbations in a coherent manner and further reduce damping. Also, the vigorous 
mixing of liquid fuel and air could generate energetic perturbations that might overcome the injector relative 
pressure-drop. Consequently, acoustic pressure oscillations can occur inside the combustor and couple with the fuel- 
air feed-system dynamics to result in non-steady heat-release that, in turn, will amplify the combustion acoustics. 
Thus, large pressure oscillations coidd occur inside the combustion chamber, disrupting flame stability and 
interfering with other subsystems (Ref. 1). Active Combustion Control (ACC) via fuel-modulation provides a 
method to suppress combustion instabilities to enable the effective operation of such lean-burning combustors for 
aero-engines (Ref. 2). 

The development of the model-based, self-tuning control for fuel-modulation to be presented in this paper is a 
part of the research on ACC at NASA Glenn to enable direct-injection lean combustion. The interested readers are 
also referred to other related reports from NASA Glenn (Refs. 3 to 5) for further information on this research, 
including a more detailed description of the research combustor (named the “NASA SNR”; Figs. 1 and 2), its 
control system, and other useful findings on combustion control dynamics with this test rig. 



Figure 1. — The single nozzle combustor. 


AIRFLOW 



Figure 2. — Schematic of the NASA SNR (in the HF configuration). 
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Figure 3. — Engine and NASA HF SNR comparison. 


I 

4 

i 

I 

< 

□ 


i 

8 


Combustor LF Instability - Run 114, Pt 55 and 42, 8/19/2000 



1 

1 

1 

1 uncontrolled 

controlled 


1 

1 

1 

i 

i 



1 

1 

1 

i 

i 



1 

1 

i 

i 



1 

1 

1 

i 

i 

i 


1 — 
1 
1 
l 

v — 

1 

1 

1 

1 

1 

1 

1 

A 1 

•\ 

\L. 

— 


0 1 — 1 i 1 — 

0 100 200 300 400 500 

Frequency, (Hz) 


Figure 4. — A case of severe “peak-splitting” under active 
combustion control for the LF instability (Ref. 13). 


In this direction, a model-based controller was developed to operate a high-frequency fuel-modulation valve 
based on dynamic pressure sensors as feedback to directly affect the instantaneous fuel-air ratio in order to suppress 
the combustion instabilities. This model-based control method for ACC (named the “Multi-Scale Extended Kalman” 
approach, abbreviated MSEK (Ref. 6)) uses multi-scale analyses (Refs. 7 and 8) and an extended Kalman fdter 
(Refs. 9 and 10) to predict and control fuel-modulation effects and the combustion thermo-acoustics. A prototype of 
the MSEK core logic (named “ the Basic MSEK ” throughout this paper) was demonstrated in 2002 on a test rig for 
an instability with high-frequency pressure oscillations (greater than 500 FIz) (Ref. 6). This NASA Glenn single 
nozzle combustion rig operates at engine pressures, temperatures, and flows using liquid jet fuel and incorporates 
many engine-like features, including and actual engine fuel nozzle and swirls, dilution cooling, and an effusion- 
cooled liner (Ref. 11). The rig (Figs. 1 and 2) was developed in partnership with Pratt & Whitney and United 
Technologies Research Center. Although the rig exhibited a higher level of background noise/disturbances than 
observed in the engine (Fig. 3), it adequately replicated an observed engine instability experience. In addition to the 
baseline configuration at “high-frequency” (F1F) the combustor rig can also be changed to a “low-frequency” (LF) 
and “mid-frequency” (MF) configuration. For the LF configuration, which exhibits ~280 FIz instability, both spool 
pieces, as shown in Figure 2 are moved to a position between the pre-diffuser and fuel injector. For the MF 
configuration, which exhibits a ~340 Hz instability, only one spool piece is moved to between the pre-diffuser and 
Riel injector. The initial investigation of the Basic MSEK with the NASA SNR in the HF configuration was 
conducted in May and September 2002. This initial investigation, performed with both the Basic MSEK and an 
adaptive controller by Kopasakis (Ref. 5), constituted the first successful demonstration of ACC on a realistic aero- 
engine combustion rig (Refs. 5, 6, and 12). Prior to that demonstration, partial success in suppressing the LF 
instability had been reported (Ref. 13). The main issue with control-performance in this and other earlier works was 
the so-called “peak-splitting” phenomenon (as illustrated by the pressure amplitude-spectra of Fig. 4). While the 
dominant instability was suppressed, other instability peaks arose around the central frequency; and consequently, 
the reduction in oscillation rms amplitude was poor. 

Such problems (perhaps due to both the valve dynamics and particular control strategies) can be overcome with a 
high-performance valve and by using variable control gains and other real-time adjustment of control parameters, as 
recently demonstrated with the “ Self-Tuning MSEK”. In this follow-on development of the MSEK, various 
parameter optimizations and “self-tuning triggers” were used to effectively minimize the time-scale averaged 
variances of the combustion sensed pressure without aggravating the existing instability at any time. Such self- 
tuning features (based on multi-scale variances or higher moments of pressure oscillations for the triggers to adjust 
the parameters and setting their increments and search ranges) are essential for handling the severe uncertainties 
inherent in combustion time-delays and other characteristics, while keeping the embedded predictor model in the 
controller simple. The resulting Self-Tuning MSEK controller showed good results on the MF and LF configurations 
of the NASA SNR, including a very fast convergence towards optimal “phase-lock.” For the LF SNR configuration, 
where the signal-to-noise ratio is high, the Self-Tuning MSEK controller was able to reduce the amplitude-spectrum 
peak by 90 percent as well as the overall pressure oscillations rms value by 80 percent without causing peak- 
splitting. 

This paper describes the improvements recently made on the Basic MSEK to create the Self Tuning MSEK. In 
sections II and III , the reduced-order model and validation results based on some open-loop forced-response data of 
the SNR are presented. The MSEK controller formulation, self-tuning features, and ACC demonstration results on 
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the SNR are presented in sections IV and V. Some conclusions and highlights of these ACC experiments, and control 
and model developments are given in section VI. 

II. Reduced-Order Thermo-Acoustics Model 

The reduced-order thermo-acoustics model for this control development, the same as presented in the previous 
report for the 2002 work (Ref. 6), is reviewed here for further clarification. Except for parameter values, this same 
model is applicable for all three configurations (HF, MF, and LF) of the SNR. It is formulated for the two dominant 
thermo-acoustics modes of the SNR as follows (see nomenclature list for the variables description): 


Xn+1 — AXn ^b u n, ^c^n, W n ) ? 

(la) 

1 0 1 
0 
11 

(lb) 

s n =K + y n 

(lc) 


In equation (la), the 4x4 state matrix A is approximated by the following expression, with A,;= exp (~2£iCO[ T): 

cos( co\T) - sin( co i T )1 r .. 

0 

sin( co \T) cos( co\T) 

cos( C 02 T) -sin( C 02 T) 
sin(w 2 T) cos( ©2 F) 

Typically, co\ (“Mode 1”) could be 271x290 radians/sec to 27ix340 radians/sec, depending on the rig; and oh 
(“Mode 2”) could be 27ix520 radians/sec or higher. The two damping coefficients C , (for i = 1, 2) are typically about 
0.01 to 0.05. The heat-release fluctuation, 0, affecting the acoustics is a 4x1 vector-valued function of delayed 
acoustics states, internal pressure feedback on fuel-flow, and fuel-modulating valve position u n (about a mean 

position); with delay effects modeled by random time-shift operators r a ; Zb , r c (of non-zero mean values) applied on 
the various variables. The calculation of 0 with these variables is described in the following intermediate steps, in 
which the parameters {/?ij} and “g” are constant, internal feedback gains; “b” represents the effectiveness of fuel- 
modulating on pressure oscillations; a reference pressure comparable to the mean pressure drop across the 
injector: 
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0 - [0 el 0 02 ]' (ig) 

The superscripted operator [...]' as in Equation (lg) and elsewhere denotes matrix/vector transpose. The imposed 
bounds on p l and p 2 following Equation (le), for the calculation of 0, was implemented differently by substituting 
the following expressions instead of 0; (for i = 1,2) into (lg): 

(2 U i l Tt) arctan (0.5 ji (0 ; + e) / U i ) (lh) 
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The inputs Q, for expression (lh) are defined as in (If), except that the bounds previously imposed on p, are no 
longer needed (Instead, the quantity 0 < [/; < I, for i = 1, 2, are constant fractions now provide the necessary upper- 
bounds for the quantity (lh)); and, “e” is a bias. These parameters as well as \h, b, and k will be chosen to match 
the open-loop forced response characteristics of the combustion chamber pressure (see section III). Note that, the 
dynamics of the heat- fluctuation term Q is very complex, as it is the result of heat convection, diffusion, air flow and 
fuel flow, and combustion. The simplistic approximation of Q in (le) to (f), with some imposed bounds (empirical, 
implemented as in (lh)), and the Euler expression for the fuel-air ratio fluctuations as in (If) is only for modeling its 
average (non-linear) effects on the acoustics. For example, the model data for the Low-Frequency Configuration as 
set-up at Glenn in 2004 are chosen as follows: 

First, choose p a = 20.0 psi; T = 0.0002 sec; K\ = 1.2 p Q 
= 24.0; k 2 = 0.4 p 0 = 8; g = 0.001; and, b = 0.3. Next, 
co\ = 1822 radians/sec; £ \ = 0.042; Ci = 1.0; h n = 1.0; 

U\ = 0.9; and, e\ = 2.64 (Mainly for matching mode-1 
characteristics). Then, ah = 3456 radians/sec; 2) = 0.01; 
c 2 = 1.0; h 2 2 = 1-0; U 2 = 0.3; e 2 = 3.2 (For mode-2 
characteristics). And lastly, h l 2 = 0.1; h 2 i = 2.0 
(The influence of mode 1 on mode 2 appeared strong 
in the LF SNR); \ = 0.0038; \ = 0.0033; 

t c = 0.0002 sec (Empirical, mean values). 

Model validation with this parameter set will be 
presented in the next section. 

“Georgia Tech valve” used for fuel-pulsing. 



ID. Model Validation With the “SNR” Test Data 

The reduced-order model and associated 
parameters for the HF configuration of the SNR had 
been validated in an earlier work (Ref. 6). 
Therefore, this section will be devoted only to 
validation of model for the LF configuration (LF) of 
the SNR. This configuration is particularly 
important for validating the model, because its 
internal instability is much more prominent 
compared to the external disturbances. The latest 
MSEK and self-tuning features were built and 
enhanced mainly with modeling data from this test 
rig configuration. The characteristics of the LF SNR 
self-sustained instabilities and open-loop responses 
to fuel-modulation by a high-frequency valve 
(Fig. 5) are shown next to validate the described 
reduced-order plant model and its parameter values. 
To highlight the agreement of model results with 
test data, this validation will be shown here with 
more refined parameters than typically done for 
preliminary control tuning prior to testing — 
preliminary tuning of control parameters prior to 
testing on the rig required only rough estimates of 
model parameters on the instability magnitude and 
the valve “ control-authority ” (i.e., the response 
characteristics of the combustor pressure under 
open-loop commands). 

Some typical non-linear response characteristics 
of the LF SNR under sinusoidal fuel-modulation 
with the high-frequency fuel-valve are shown in 
Figures 6 and 7. In these figures, the signal named 
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Figure 6. — (Upper subplots) The LF SNR self-sustained 
instabilities (oscillations in psi); (Lower) Effects 
of 100 Hz sinusoidal valve commands with 3-V 
amplitude on the instabilities, resulting in about 
10 psi change in pressure oscillations. 
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Figure 7. — Effects of 200 Hz sinusoidal valve commands of 3-V 
amplitude on the LF SNR instabilities, resulting in about 
7 psi change in pressure oscillations. 
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Simulation of SNR LF, responses to 100Hz, 3-Volts valve excitations 



Figure 8. — (Upper subplots) Model results of LF self-sustained instabilities; (Lower) Effects of 100 Hz sinusoidal 
valve commands of 3-V amplitude on the LF SNR Model instabilities. 


“PLA1C1” is the combustion pressure oscillation signals sensed at about 2 inches downstream of the fuel injector; 
and “REFAread” (in Volts) represents the valve command. For example, the two upper subplots of Figure 6 show 
the amplitude spectra and time -do main trace of the combustor pressure self-sustained oscillations (in psi) without 
fuel-modulation. The main instability frequency is about 315 Hz. The higher frequency mode, at about 620 Hz, is 
also visible on the amplitude spectra plot. The two lower subplots of Figure 6 show the responses of the internal 
pressures to 100 Hz, 3-V sinusoidal open-loop commands on the valve. The lower subplots in Figure 6 as well as 
Figure 7 (for 200 Hz open-loop inputs) clearly show the nonlinear response characteristics of the combustion 
instabilities and the significant effects of the fuel-valve pulsation. 
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Figure 8 shows the simulation data of the LF SNR 
reduced-order model to be compared to the test data of 
Figure 6. It shows that, the modeled self-sustained 
instabilities as well as their responses to the 100 Hz, 3-V 
sinusoidal fuel-modulation open-loop commands, 
including non-linear characteristics, compare well on both 
the frequency-domain and the time-domain with the test 
data of the rig. Likewise, the model results of Figure 9 
compare well with the test results shown in Figure 7. 

IV. The Self-Tuning MSEK Controller 

The top-level structure of the Self-Tuning MSEK 
controller for ACC is still the same as that of the “basic 
MSEK”(Ref. 6) tested in 2002 on the HF SNR (see 
Fig. 10). The control scheme formulation, and additional 
self-tuning features as in the latest version will be 
described in more detail in the following sub-sections. The 
subsection titles refer to the individual sub-block names in Figure 10. Except for the last subsection ( Automatic 
Timing of Critical Control Parameters), these subsections are reviews of the basic control scheme that had been 
briefly described in the previous work (Ref. 6). 

A. Multi-Scale Tone Analysis Using “QA” 

This pre-filtering part of controller MSEK (Fig. 10) is performed by a “Quadratures Filter Bank” of wavelet-like 
filters to extract the dominant “tones” of the sensed pressures. By definition, a “ Quadrature Analysis ” (QA) on a 
frequency range of a signal stream {s„} (a time-series with n being the time index) is any decomposition of the 
signal contents within that range into two mutually orthogonal components {Si,,} and { S 2n } by using an 
Orthogonal Wavelets Decomposition and Partial Reconstruction scheme (e.g., using “quadrature mirror filters”) 
(Ref. 7). For simplicity, instead of such a wavelet scheme, a pair of Finite Impulse Response (FIR) band-pass filters 
is designed for each frequency range to extract (approximate) “quadratures” from the signal: One filter with unity 
gain and 45-degrees lead at the center of the frequency range of interest, to extract the time-series { S\ „ }; and 
another, also with unity gain, but with 45-degrees lag at the same frequency, to extract { S 2n }• Using the described 
QA process, for each of the two main frequency ranges (LF and HF) of the SNR, a “filter-bank” consist of three 
Quadrature Filter pairs (hence a minimum of six filters for each range) was used to analyze the sensed pressure 
oscillations and, for each time point /?, extract the dominant pair among the three pairs of quadrature components. In 
particular, for the LF range, one pair of Quadrature Filter as described is designed for the central LF, and the other 
two are for 0.98 LF and 1.02 LF, respectively. These are then used to extract the dominant LF QA sequence, 

denoted { S „ = [5 1 „ S 2 „ ] ' } as in the formulation of the predictor dynamics in the next section. Likewise is done 
to extract the HF QA sequence which will be denoted { S „ = [S S 2n ] }. The advantage of this (approximate) 

QA process is its simplicity; also, the polar angle of the vector [ S\„ S 2n ]' (being the phase angle of { S\ n } ) 
represents the instantaneous “phase” angle of the tone i plus n/4 radians. 

B. EK States Prediction 

i 

Although the evolution equation for S „ can be derived from the state-equations (la) to (lh) to develop the 
prediction scheme, such route is avoided because it is excessively complicated for real-time calculation. Thus, the 
state-equations (la) to (lh) will be used directly in the following development of the Predictor (Eqs. (3a) to (3e)). 
Note also that, if the “QA” process is applied to the output {s„} simulated by the model in equation (lc), then the 
resulting S „ is approximately C- X‘ n (for i = 1, 2) with about 45-degree phase-lead. Hence, the following fictitious 
random process would be a good projection of the actual state variables X„ ahead N sampling steps (such that, N x T 

i 

~ total delay-time) if certain augmented input “L can be designed based on S „ (for i = 1, 2) to make the output 
y n track the measurable output “ c x X\ + c 2 X}, ” despite the random disturbances w and v: 


Simulation; SNR LF Config, responses to 200Hz, "3-Volts" valve excitations 



Frequency (Hz) 

Figure 9. — Response characteristics of the 
LF SNR Model to 200 Hz valve 
commands of 3-V amplitude. 
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Figure 10. — The MSEK controller functional block diagram. 
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The random effects w and v need to be included here in the design of the Kalman filter for q n estimation, so that 
the predictions would be robust against disturbances. Also, to simplify the Kalman filter, the projected “g”- term 
shown in Equation (2a) will be approximated using S „ (for i = 1, 2) (but, without time-delays in the formulation; 
(see Eqs. (3c) and (3d)) rather than by integration as part of the above process. Furthermore, the first two terms of 
the right-hand-side of Equation (2b) add up to Mq n _ N , where M = [c] /, q/] , and “7” denotes the 2x2 identity 
matrix. So, if the 5th-order Pade approximation is used for the N-step delay to apply on the 2-vector M q n to 
compute the delayed sequence, the combined state matrices of that transformation individually applied on the two 
vector components would result in the state-space matrices U, F, E, D (of dimensions 10 x 10, 10 x 2, 2x l(), and 2x2, 
respectively) that yield the following reformulation of Eq. (2b): 


n + 1 = U z n + FM q» 

(2c) 

y n « Ez n + DMq n + v 

(2d) 


Equations (2a), (2c), and (2d) result in the following Quadratures Kalman Predictor obtained by estimating q n ; 
the estimation is denoted q n for simplicity, and meant as a projection for X H , N time-steps ahead: 


= Aq n + Q + K x (e Z n + DMq „ - Sj - S„ 2 ) 


9 n+ 1 

Z n+ 1 = UZ n + FMq n + K 2 (eZ„ + DMq 


s») 


(3a) 

(3b) 


Pi 

Pi 


k u 0 k l2 0 


v 2 1 


0 k 22 0 


y 


subjected to : p\ mn < /f , p 2 < p\ 


(3c) 
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Qi = K i 



1 , 2 ) 


(3d) 


Q = [ 0 Q\ 0 fi 2 ]' 


(3e) 


In Equation (3c) the parameters “ k\ j ” are related to “ h\ j ” such that: /?[ j ~ c, j (for i,j = 1,2), where “c” and 
“ /7i j” are the model parameters shown in Equations (lb) and (le). The matrices K\ and K 2 (of dimension 4x2 and 
10><2, respectively) are the two vertical sub-blocks of the Kalman matrix computed for states estimation, assuming 
certain covariance matrices for the random vectors v and w. Note that, the time-delay of the “Q”- term in the actual 

process X n was dropped to affect the desired lead time for this prediction scheme (i.e., comparing Eq. (3c) to (le)). 
It is also sufficient to use the simple bounds after Equation (3c) instead of the more refined expressions (lh). 

C. The Damper 

If S = { [.S' i „ , S 2 „] f } denotes the time-series of “quadrature pairs” extracted from a signals { s„ } by a “QA” 
process as described previously, then the linear combination (associated with an angle a, in radian) denoted 
S a = {-S' i n cosa - S 2 n sina} represents the extracted tone of {s„} phase-shifted by “a+jr/4” radians. Such 
combination is used on the predicted quadratures, namely S l n = c ; - q,\ (for i = 1, 2), to prescribe a Linear Quadratic 

Regulator-type damper (This part of the fuel-modulation command, as in Fig. 10, is distinguished here by the 
superscript “d,” with the time-index suppressed and the superscript being the mode index): 


= d 


l S\ 


+ M )^ 2 


(4) 


where is a chosen “ mode-1 fraction ” (with 0 < / < 1) (basically a Loop-Shaping means for performance 
optimization), and “p” is the overall gain for the “Damper”. The terms SL (for i = 1, 2) are the phase shifted signals 

Pi 

constructed from the predicted “quadratures,” using the respective phase angles //, = fi + e)j n (for i = 1, 2), with [i 
being a certain bias and e>; „ , a predicted phase-drift. The drifts Si n (over the time interval of expected control 
transport-delay, for tones 7=1,2) are predicted by a dynamic, linear regression scheme applied on the quadrature 
phases { ai n } as defined in equation (5) of the next sub-section. 


D. The “Phase-Adjusted Reconstruction” and Tone Suppression 

This part of the control scheme is for countering the non-steady heat-release effects. For this purpose, the 
instantaneous amplitudes and phases of the combustion pressure near the LF and HF are first estimated by projecting 

the quadratures pair Sj, (for i = 1, 2, as 2x1 vectors) onto an orthogonal frame rotating at the respective speed cp 
(referring to the frequencies in Equation (Id)). To minimize the effect of random noises, these projections are next 
passed through low-pass filters with cut-off frequencies at 100 and 200 Hz, respectively, before computing the 
amplitudes and phases. Note that, these cut-off frequencies are about 1/3 of the respective central frequencies (So 
that variations in amplitude for every three oscillation at such frequency should be detected). Those estimations form 
the following expression for feedback to the tone-suppression logic (for modes 7 = 1,2): 


a : = r:,, cos 


(« 


CO I 


T + 


(5) 


where { /•, „ | is the computed amplitude, and { a, n } the phase angle at time t„ of the filtered, projected quadratures 
vector. Then, { a x „ | is used to predict the phase drifts S t n (for the expected control transport-delay, for 7=1,2) 
using a dynamic linear regression scheme. Since the quantities cr, should be highly coherent with the non-steady 
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heat-release, the fuel-modulation command u n will include a compensation term (distinguished by the superscript 
V) of the following form for Tone Suppression (see Fig. 10) to minimize these effects: 


= y 


hln 


1 n 


cos \n co 


T + 


'■ 1 n 


+ a 


+ 8 


In 


-( 1-0 


X 2n r 2 n cos \ n co 2 T + a 2n + a 


+ 8 


2 n 


, (6) 


where “y” is the overall gain for Tones-Suppression and “0 < / < 1”, the same “ mode-1 fraction’’ used in 
Equation (4). A varying factor, { y in } (for i = 1 , 2), is also introduced to allow for gradual reduction of the Tone- 
Suppression actions when “ r, ” is sufficiently small (This is essential for rendering absolute stability to the 
controlled model as shown in the 2002 study) (Ref. 6). It is basically a unit-gain with an integral adjustment, 
calculated by digitizing the following expression (also, n ” is subjected to certain bounds y min and y max to 
maintain numerical stability): 


Xi= 1 + h | { r i ~ pM , ( 7 ) 

where “ tj ” is a constant, and “p-, ” is a reduction level expected to be maintainable for “ r\ ” with this scheme. In 
summary, the Tones-Suppression command is generated with these gains on “ rr, ” after having phase-shifted it by a 
bias d plus §\ n . 

E. Automatic Tuning of Critical Control Parameters 

The fuel-modulation command { « n } issued to the valve is the sum of the expressions (4) and (6). For all SNR 
configurations, the parameters in the predictor as well as in the described command scheme appear to be fairly 
deterministic, except for the phase biases ( a and ft )and the control authority represented by “ b ” (see Eq. (3c)). The 
critical control parameters of the MSEK that require automatic tuning are: a , /? , b , p , and y. The remaining 
parameters, including (also of importance) can be preset based on the model, or calibrated empirically (In 
particular, the parameter can also be searched in real-time automatically, at a very safe and slow pace as 
demonstrated on the model as in the last subsection of section V (Control Evaluation). 

The three parameters a , ft , and b are determined by the Parameters Tuning logic (Fig. 10) via sequential, global- 
type searches. These searches are to minimize the Time-Scale Averaged Pressure Variance (in Fig. 10, obtained by 
applying a 0.2-Hz low-pass filter on the total square-norm of the two vectors S l „ , for i = 1, 2). The searches are 
activated only if performance degrades below a certain threshold (a scalable threshold based on the model). The 
search neighborhood is also narrowed down if performance is improving. In parallel, fast management of p, and y 
(with a common variable factor) are done to avoid aggravating the instabilities at wrong values of a , P , and b . 
That is, these gains are dropped quickly to zero when an overshoot trend is detected in the pressure variances near 
the instability frequencies. 

The “ Basic MSEK ” (the core of the “Self-Tuning MSEK,” namely the control actions formulated in 
Equations (4) and (6) was implemented for testing in 2002 on the HF SNR, but with only manual activation of the 
real-time search for optimal a , f , and b . Also, the gains p and y were probed and adjusted manually during the 
tests in 2002. But, for the Self-Tuning version of MSEK, multi-scale based variances and third-moments of pressure 
oscillations are used as activation functions in the Parameter Tuning block to automatically trigger the searches and 
fine-tunings as described. Upon search activation, the search diameter is adjusted according to a Time-Scale 
Averaged Pressure Variance computed with a sufficiently large time-constant for averaging. Also, the gains “//” and 
“y” in the Self-Tuning version are adjusted automatically, and continuously based on a Time-Scale Averaged 
Pressure Variance computed with a much smaller time-constant for averaging. 

V. Controller Evaluation 

First, the “ Basic MSEK ’ (the core logic of the “Self-Tuning MSEK”) was validated and successfully 
demonstrated at UTRC in 2002 with the HF SNR, also using the high-frequency valve designed by the Georgia 
Institute of Technology for fuel pulsing (Fig. 5) (Ref. 6). The validation for this part of the controller will be briefly 
summarized in the next subsection. 
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A. Model-Based Validation and Experimental Results of the “Basic MSEK” 

The SNR was operated at a pressure, temperature, and fuel-air ratio corresponding to a mid-power engine 
condition (T3 = 770 °F, P3 = 200 psia, fuel-air ratio = 0.03). This HF configuration of the SNR (Fig. 2) exhibited 
roughly the same engine combustor instability behavior as shown in Figure 3. The SNR combustor pressure was 
sensed about 2 inches downstream of the fuel injector. The control algorithms were implemented on a real-time 
processor for an I/O rate of 5000 Hz. As in the LF case, the parameters of the (reduced-order) model for that HF 
configuration (to be shown in the next paragraph) are with “psi” as unit for the pressure feedback, and Volt for the 
fuel-modulation command u a issued to the fuel-valve. This valve command is limited to 2 Volts in amplitude for the 
HF rig, which corresponds to the maximum valve displacement allowed about the mean positions. The mean valve 
position was set for the optimal flow number regarding control-authority. 

Without loss of generality, the measurement gain and the internal feedback gain for the dominant mode (Eq. (lb) 
and (le), respectively) were normalized and balanced such that: C 2 = /i 2 2 = 1-0. With this balance normalization and 
the mentioned unit choice, the rest of model parameters were then chosen to match empirical open-loop 
characteristics of the combustor (e.g., Fig. 3). And, this was done in the order of matching mode 2 characteristics 
and responses to fuel modulations first, as it is more dominant, and then mode 1 : 

First, p 0 = 20.0 psi; T = 0.0002 sec; K\ = k 2 = 0.6 p 0 = 12.0; g = 0.01; and, b = 0.2322. Next, ah = 
3456 radians/sec; C, 2 = 0.04; and, p’ mSK = 0.09 p 0 (for matching mode 2 characteristics). Then, oj\ = 1822 radians/sec; 

= 0.3; Ci = 0.8; h u = 1.0 (for mode 1 characteristics); and h\ 2 = h 2 1 = 0.1 (Modal-couplings in the HF 
configuration of the SNR apparently were low); x a = 0.0038; f b = 0.0033; x c = 0.0002 seconds (Empirical, mean 
values). 

Figure 1 1 shows the predicted Amplitude Spectra and Probability Density of combustor pressure oscillations; the 
control-ON compared to the control-OFF case for a rather high level of disturbances comparable to the actual HF 
SNR environment. Figure 12 shows similar results from test data in 2002 with the Basic MSEK controller on the HF 
SNR. Both the model prediction and test data show a reduction of about 30 percent in pressure oscillation amplitude 
at the HF peak. Reference 6 also includes some analytical results showing the controlled “absolute stability” 
achieved by this core logic. Note that, controlled “absolute stability” is possible if the time-delay of fuel-modulation 
effect (“t 4 ”, Eq. (la)) are comparable to, if not less than, that of the internal heat-release feedback (“x a ”). 


Amplitude Spectrum - MSEK ACC simulation at high random disturbances 




Closed-loop 



/v\ 

Open-loop 










°-4 - 3 - 2-101234 

Pressure in PSI (Probability Density Distribution) 


Figure 11. — MSEK, HF SNR Simulation: pressure 
oscillations amplitude-spectrum (Upper) and 
probability density of pressure oscillations 
(Lower), with/without control. 



Figure 12. — (9-1 1-02 HF SNR test) Amplitude spectra 
showing the effects of the Basic MSEK controller on 
the combustion instability peak pressure oscillation. 


NASA/TM— 2006-213855 


11 


B. Model-Based Validation and Experimental Results of the “Self-Tuning MSEK” 

The “Mid-Frequency” (MF) and “Low-Frequency” (LF) SNR configuration were also used in the 2004 ACC 
demonstrations at NASA Glenn. The MSEK and it self-tuning features and parameters were designed and validated 
on the Reduced-Order model shown in sections II and III , and also on a Sectored 1—D model (1-D CFD, by Paxson, 
Ref. 14) before testing on the SNR (see Figs. 13 and 14 for some validation results of this controller on the Sectored 
1-D model). 

The main findings for the optimal operation of the Self-Tuning MSEK on the LF rig include: The setting of the 

valve “flow-number” for optimal fuel-modulation effects; the nominal search steps (for a , f , and b ); the threshold 
for initiating global search; and lastly the optimal loop-shaping fraction “/”. Initially, “/” had been set at 0.5. With 
that value of and the optimal search steps and threshold, the automatic tuning actions of the MSEK controller 
was able to suppress the instabilities by 80 percent after having been engaged for about 5 sec. Flowever, such 
suppression could not be held steady with setting / = 0.5 because of some low-frequency oscillations (see Fig. 15). 
Such low-frequency oscillations are most likely due to the mutual “ beating ” between mode-1 instability (315 Flz) 
with the fuel-modulation effects by the controller near 320 Flz that appeared to have suppressed mode-1 almost 
completely. The left subplot of Figure 15 shows the amplitude-spectrum of the pressure oscillations of the middle 
subplot for this case; and the right subplot shows the spectrum of the valve command signal. This command 
spectrum is broadband, but concentrated mostly around 300 Flz. This setting reduced the peak of the pressure 
amplitude spectrum down to about 1.2 psi/Flz. This is a significant reduction of the instability, as compared to when 
the controller is turned off, which is about 1 1 psi/Flz (e.g., the case of Fig. 6). 

When is set near 1.0 such as in the case of Figure 16, the resulting command spectrum was sharper and more 
concentrated at 300 Flz. But, the instabilities appeared to be suppressed by only 70 percent and yet not steady, again 
because of some low-frequency oscillations. Although the original instability that had been at 315 Hz was 


Amplitude Spectrum - ACC slm 



Figure 13. — MSEK on a HF Sectored 1-D model. 



Figure 14. — MSEK on an LF Sectored 1-D model. 


Run 420, 040519 , Point pt097 Run 425 > 040602 , Point pt067 




Figure 15. — Preliminary test results of the self-tuning 
MSEK on the LF (315 Hz) SNR-May 19, 2004. 


Figure 16. — Preliminary test results of the self-tuning 
MSEK on the LF (3 15 Hz) SNR-June 2, 2004. 
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eliminated entirely this time, the fuel-modulation produced an undesirable disturbance at 273 Hz with unsteady 
amplitude. 

The instability suppression by the Self-Tuning MSEK was best and steady when the command was more 
broadband towards the high-frequency range as shown in Figure 17. This was when was set near 0.1. Under this 
setting, the controller was able to suppress the oscillation (rms) amplitude from about 10 psi (e.g., shown in the 
upper subplots of Fig. 6) to about 1.7 psi. This was so because the instabilities were concentrated mostly at 315 Hz 
while fuel-modulation with the GT valve was able to affect this frequency most deeply over a fairly wide range of 
commanded frequencies. Apparently, fuel-modulation at high-frequencies for the LF configuration would not excite 
the higher modes, but instead affects significant energy transfer towards suppressing the low-frequency instability. 
So, the role of “ 0 < / < 1 ” is clearly that of typical loop-shaping with a sensitive control parameter. This seems to 
be essentially the same phenomenon discovered previously by Kopasakis on the same rig when focusing control on 
the “2nd harmonic ” with an adaptive phase-shifting controller (Ref. 3). 

Lastly, the MF SNR exhibited an instability at 340 Hz (Fig. 18(a)), a frequency fairly close to the LF; but, its 
nature is more like that of the HF instability. The plots shown in Figure 18(b) are from a preliminary study of the 
MSEK on the MF configuration that show some moderate reduction in instabilities; about 50 percent reduction of 
the amplitude spectra peak relative to the noise floor compared to open-loop. Also, the rms amplitude of overall 
pressure oscillation was reduced to about 1.4 psi. This is probably the limit for active fuel-modulation effect with 
this valve (Fig. 5) for this case. This is consistent 
with the 2002 results on the I IF rig, and the test 
results on the LF rig shown in Figure 17 (Noting 
that, the time-plot in Fig. 17 also showed peaks near 
4 psi, comparable to the time-plot of Fig. 18(b)). 

C. MSEK Model-Based Control Dynamics 
Discussion 

The LF Configuration experimental result of 
Figure 17 is an important case that had not 
previously been simulated. It was found from 
simulating this case that the non-linear heat-release 
over the frequency range between mode-1 and 
mode-2 seems to be crucial for the steady 
suppression of the instability. Therefore, some 
numerical results of the MSEK on the validated 

Figure 17. — The LF (315 Hz) suppressed by the self- 
Tuning MSEK controller-June 2, 2004 Test. 


Run 425, 040602 , Point pt080 



PLAICIpsi 


Run 429, 040707 , Point pt041 




Time, sec 



Figure 18(a). — MF self-sustained instability. 


Figure 18(b). — MF instability suppressed by MSEK. 
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reduced-order model are presented here for further 
understanding of the control dynamics and for validating 
the model against test data for the close-loop case. First, 
Figure 19(a) shows the self-tuning control in action on the 
LF model, and the combustion pressure being quickly 
reduced to a level comparable to an actual test case on the 
LF SNR (Compare Fig. 17 of an actual test case with a 
“zoomed-in” portion of the simulation run case in 
Fig. 19(a), as shown in Fig. 19(b)). The “ mode-1 fraction, 
F of Equations (4) and (6) of the Damper and Tone 
Suppressor actions, used as a “ loop-shaping ” parameter, 
is shown in the third subplot of Figure “19(a)”. Note that, 
comparing to those results within the first second time 
interval when the control had not been activated, the 
pressure oscillations afterwards have never been 
worsened by the tuning actions. Figure 19(c) shows the 
difference in the unsteadiness of heat-addition of two 
zoomed-in time portions in Figure 19(a): one portion at 
0.5 sec from start when the controller had not been 
activated; and another at time 20 sec when the controller 
parameters has already converged to a good region. It 
shows that, fuel-modulation under good control likely has 
the effects of preventing the heat addition from dropping 
too low in resonance with the acoustics. Thus, fuel- 
modulation might also be effective for stabilizing the 
flame when the combustor is operated near its “ lean - 
blowout” margin. Such potential application of active 
fuel-modulation (e.g., to “lean-direct injection 
combustors), however, needs further investigations. 

VI. Conclusion 

This combustion stability control scheme, the self- 
tuning Multi-Scale Extended Kalman (MSEK) controller, 
has been used to experimentally suppress thermo-acoustic 
instabilities up to 90 percent. The controller employed a 
Quadratures Filter Bank to perform multi-scale analyses 
in real-time on the sensed pressure signals for control 
feedback. An approximate evolution equation of the 
extracted “Quadratures” (representing the acoustics 
states) was used to formulate a Kalman predictor for these 
acoustics states for the damper logic. The variable gains, 
phase bias, and other automatic control tuning logics were 
based on such multi-scale analysis. This model-based, 
self-tuning control scheme was thoroughly tested on a 
liquid-fuel combustor rig operated at engine-like 
conditions, for three configurations emulating different 
cases of combustion instability. The combustion control 
study on this rig confirmed earlier NASA observations 
and provided a plausible explanation about an important 
internal mechanism of cross-frequency couplings in 
combustion thermo-acoustics. The dynamic couplings are 
very prominent in its Low-Frequency configuration, and 
in general could be strongly affected by local fuel/air 
modulations. The reduced-order thermo-acoustics model 
for control development, using a simplistic parametric 


A Simulation of the Self-Tuning MSEK on the Reduced-Order Model 




Figure 19(a). — Simulation results with the self- 
tuning MSEK controller on the reduced-order plant 


Simulation: Pressure Spectrum Pressure Oscillations Command Spectrum 



Frequency (Hz) Time (Seconds) Frequency (Hz) 

Figure 19(b). — Simulation analysis of the self- 
tuning MSEK controller (see Fig 19(a)). 


Simulation: Pressure Oscillations Pressure Oscillations 




Heat Unsteadiness Frofile 



■ 3 1 ■ 1 
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Figure 19(c). — Instability and heat unsteadiness — 
simulation (see Fig 19(a)). 
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formulation for the non-linear heat-fluctuations that drive the acoustics, exhibits the same response characteristics of 
the rig on both the frequency-domain and time-domain. The self-tuning MSEK based on this model was able to 
adjust it own parameters and control gains to suppress nearly completely the combustion instabilities by modulating 
the fuel injector flow with a high-frequency valve. The main contributors to the effective suppression of combustion 
instabilities as demonstrated by this controller were the control-authority of the high-bandwidth fuel valve and the 
cross-frequency couplings in heat-release under fuel modulation. 
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